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Using the Continuous Time Quantum Monte Carlo Loop algorithm, we calculate the temperature 
dependence of the uniform susceptibility, and the specific heat of a spin-1/2 chain with random 
antiferromagnetic and ferromagnetic couplings, down to very low temperatures. Our data show 
a consistent scaling behavior in both quantities and support strongly the conjecture drawn from 
the approximative real-space renormalization group treatment. A statistical analysis scheme is 
developed which will be useful for the search scaling behavior in numerical and experimental data 
of random spin chains. 

75.10.Jm, 75.10.Hk, 75.10.Nr 



Spin chains have for many decades attracted much in- 
terest in condensed matter physics and statistical me- 
chanics. In recent years random spin chains have been 
investigated intensively because they provide a simple 
model to study the interplay between quantum effects 
and disorder. There are various realizations of quasi one- 
dimensional random spin systems in nature. One class 
is represented by compounds like Sr^CuPti^xIrxOg 
y. While the pure compounds Sr^CuPtOe (x = 
0) and Sr^CuIrOe {x — 1) are antiferromagnetic 
(AF) and ferromagnetic (FM), respectively, the alloy 
Sr^CuPti-xIfxOf, contains both AF and FM couplings 
whose fraction is simply related to the concentration x 
of Ir. A corresponding minimal model is the spin-1/2 
Heisenberg chain, where the nearest neighbor exchange 
coupling is + J or —J with certain probabilities [0. An- 
other related example is realized in the low-temperature 
regime of randomly depleted AF spin-1/2 Heisenberg lad- 
ders. Two-leg Heisenberg ladders have a resonating va- 
lence bond ground state which has short-range singlet 
correlation and a spin excitation gap [|| . If spins are de- 
pleted at random the low-energy properties are changed 
drastically. Each site which lost a spin is accompanied by 
an effective spin 1/2, and the residual interaction among 
these spins is randomly FM or AF with a wide distribu- 
tion of coupling strengths [Q,|| . 

In this letter we focus on random FM-AF spin chains 
with the full spin rotation symmetry. One of the most 
powerful techniques to study such systems is the real 
space renormalization group (RSRG) method [p|-fi2[. 
This method was introduced to study random AF spin- 
1/2 chains [p[-pi)[, and was recently adapted to the study 
of the class of systems introduced above [ pTp^ . The ba- 
sic idea of the RSRG method is the iterative decimation 
of degrees of freedom by integrating out successively the 
strongest bonds in the spin chain. In this way the distri- 
bution of coupling strengths and spin sizes is renormal- 
ized. For many cases a universal fixed point distribution 



is reached in the low-energy limit. The random AF spin- 
1/2 chain belongs to a universality class different from 
that of the random FM-AF spin chain. The former is 
characterized by a random singlet phase where each spin 
tends to form a singlet with one other spin even over 
very large distances. For the FM-AF spin chain, how- 
ever, spins correlate to form effective spins whose average 
size grows with lowering of the energy scale |ll[-p^ . A 
result of the RSRG method is that the average number 
of original spins included in a single effective spin for a 
given (energy or) temperature T scales as I ex T^'^" , and 
the average spin size S ex T"" for T ^ [0[2). In the 
latter case this RSRG treatment gives a consistent de- 
scription of the random spin chain, but relies on certain 
approximations which have not so far been independently 
tested. 

Here we use a different approach to investigate the 
random FM-AF spin-1/2 chain and its scaling behav- 
ior. Applying the recently-developed Continuous Time 
Quantum Monte Carlo (CT QMC) loop algorithm 0, 
we calculate the temperature dependence of the uniform 
susceptibility and the specific heat. The powerful CT 
QMC algorithm enables us to simulate systems down to 
very low temperatures, and thus to observe the signs of 
the low-energy scaling regime of the fixed point. We 
find excellent agreement between our analysis and the 
RSRG treatment by Westerberg et al. ||ll|,|l|l. Beyond 
the consistent evaluation of the scaling exponent a, our 
simulation gives various quantities over a wide range of 
temperatures and our discussion provides a technique to 
analyze (numerical or experimental) data for the low- 
temperature scaling regime. 

For the numerical simulation we use as the model of a 
random FM-AF spin-1/2 chain 
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where Si is the r spin of the chain and Ji, the cou- 



pling strength of the i*'' bond, takes random values of 
both positive and negative sign. We assume the bond 
distribution 
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where Jq is the maximal coupling setting the energy scale. 

According to Westerberg et al. [|l2|, the low energy 
behavior of the random FM-AF chain is independent of 
the initial distribution as long as jy^P{J) is regular for 
J — > with He ~ 0.7. Therefore, assuming a particu- 
lar distribution (Eq. (g)) is not a restriction and our 
results apply qualitatively for all random spin chains of 
this class. 

Using the CT QMC method, we simulate 100 samples 
of spin chains of 400 sites in a temperature range down 
to temperatures as low as T = Jq/IOOO. In the CT QMC 
loop algorithm the calculations are performed directly 
with a Trotter time interval At — 0, so no extrapola- 
tion in the Trotter time interval is necessary. To obtain 
good statistics, we consider 100 different samples of the 
spin chain following the distribution in Eq. (^. The 
temperature dependences of the physical properties are 
calculated for each chain separately, and then averaged 
over the 100 samples. 

An important result of the RSRG approach is that 
the uniform susceptibility approaches the Curie-like be- 
havior in the limit T ~* 0. It is possible to determine 
the limiting Curie constant from the initial spin size dis- 
tribution and the ratio of FM and AF bonds. The low- 
temperature deviations from this Curie constant give im- 
portant information on the scaling properties mentioned 
above. For the scaling analysis we calculate the prod- 
uct Tx{T) with the QMC simulation over a temperature 
range, < T < Jg as shown in Fig. 1. According to the 
RSRG scheme with lowering temperature gradually the 
spins correlate within clusters of growing length I. Each 
such cluster represents an effective spin. This leads in 
our case to an effective diminishing of the spin degrees of 
freedom such that Tx{T), as a measure for the still un- 
correlatcd spins at a given temperature, decreases mono- 
tonically from the large-T value of 1/4, the free spin-1/2 
limit. The linear dependence in the intermediate tem- 
perature range reflects the initial bond distribution P{J) 
in Eq. (0) . This can be very clearly observed in the com- 
parison of the QMC results with the finite size systems, 
L = 2, 3 and 4, for which we have performed the exact 
bond average (Fig. 1). The deviation from this linear 
behavior for T < O.OSJq indicates the turn towards the 
low-temperature scaling regime. Here the effect of the 
renormalization of the distribution function becomes vis- 
ible. Note that such renormalization effects are already 
possible for finite size systems with L = 4 as can be seen 
in Fig.l. In the low-temperature regime thermodynamic 
quantities should approach universal behavior which we 



will now analyze based on the scaling assumption of the 
RSRG method. 
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FIG. 1. Temperature dependence of the uniform Curie con- 
stant of random chains of different length L with a bond dis- 
tribution according to (Q). The data for chains of two, three 
and four sites are calculated exactly, integrating over the dis- 
tribution (|) for each bond. The data for L = 400 are QMC 
results. 

The fact that the correlation of spins occurs within 
definite clusters permits the use of a special statistical 
approach for the discussion of the thermodynamic prop- 
erties. The effective spin of a cluster of length I is equal 
to the ground state spin quantum number. For spin-1/2 
degrees of freedom coupled randomly by FM or AF bonds 
the probability distribution for these quantum numbers 
is given by. 
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where S ranges from to 1/2 {S: integer for I even 
and half- integer for I odd). Because the system contains 
no frustrating couplings, the spin quantum number S is 
identical to the total spin of the corresponding completely 
correlated cluster of classical spins [^,|ll|,|l2[ . 

The average value of the Curie constant per site for a 
finite, completely correlated cluster of length / is 
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Now let us consider the distribution of cluster lengths I 
for a given temperature. It is natural to assume that the 
bonds freeze uncorrelated among each other, i.e. each 
bond is frozen with a certain probability independent of 
the location within the chain. Thus the distribution of 
I, Pi{l) has an exponential form. As mentioned above, 
the RSRG approach predicts a scaling behavior for the 
temperature dependence of the average cluster length size 
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where A is a dimensionless proportionality factor. This 
can now be used to calculate the low-temperature be- 
havior of Tx- If ^ ^ 1, it is justified that we use the 
distribution functions pi{S) and Pi{l) in their continuum 
approximations 
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for Z > 1. (We have compared the following results also 
numerically with the discrete forms and found only minor 
deviations even for small /.) The averaging over both 
distributions leads to 
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1 [-1/2 



12 6V2 
1 T" 



0{l 



-3/2N 



12 6Ax/2 



+ 0{T''^), 



(7) 



where for the third line I was substituted by Eq. (|^). 
Here the idea enters that the effective spins formed by 
each cluster behave independently under thermal fluctu- 
ations. 

Fitting the leading two terms of the above form (Eq. 
(M)) to the QMC data gives an estimate of the scaling 
exponent a and the proportionality factor A. Fig. 2 
shows Tx as a function of T" for different a. Using the 
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FIG. 2. Uniform Curie constant per spin as a function of 
T" for different a. The best linear behavior (Eq. (W)) is found 
for a = 0.21 ± 0.02. The sohd line shows a fit of form (m) for 
a — 0.21 and fixed intercept 1/12. This leads to a value of 
A = 1.0 ±0.1 . 
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FIG. 3. Specific heat divided by temperature T as a func- 
tion of T. The circles are the QMC results, while the solid 
line was calculated as the derivative of Eq. (hd) with respect 
to T. The inset shows the QMC results for the specific heat 
as a function of temperature over a large temperature range. 
The error bars are smaller than the symbols. 

knowledge of the exact zero-temperature value of Tx 
(= 1/12) we obtain from this plot a = 0.21 ± 0.02, in 
good agreement with the RSRG result (0.22 ± 0.01) [0. 
The proportionality factor is determined as A = l.OibO.l. 
These results confirm that our QMC results for chains 
of length L = 400 are not affected by finite-size effects, 
because according to Eq. (0) the average length of the 
correlated clusters, even at the lowest simulated temper- 
ature Jo /1 000, is only Z ss 18, and therefore much smaller 
than L. 

Now we turn to the specific heat C, which is deter- 
mined from the numerical derivative of the internal en- 
ergy calculated in the QMC simulations. The errors in 
the specific heat increase significantly with lowering tem- 
perature. Nevertheless, C could be determined reliably 
down to T « Jo/150. The QMC results for the specific 
heat and the ratio C'/T as a function of T are shown in 
Fig. 3. 

Since C'/T is the derivative of the entropy a with re- 
spect to T, we may compare it to a calculation of a by 
our statistical cluster analysis. Again we assume that for 
given temperature each correlated cluster behaves inde- 
pendently. Hence a cluster of length I contributes 
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to the entropy. Here it is necessary to use the discrete 
expressions for the distribution functions because large 
deviations arise in the continuum approach using Eq. m). 
The discrete distribution pi{l) now has the form 
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Note that this distribution is normahzed for a sum over 
I from 1 to cxD. The entropy per site of the infinite chain 
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where I is given by Eq. (pf). Using the values of a and 
A determined above, the entropy and hence C/T may be 
calculated from Eq. ( [lO| ) . The result is shown as the solid 
line in Fig. 3. We find good agreement with the QMC 
results for the low-temperature range and emphasize that 
no further fitting was required to obtain this result. At 
higher temperatures this approach must fail, because the 
scaling behavior of I (Eq. (|^)) is no longer valid. 

As a further test of the consistency of our analysis, we 
estimate the area below the curve C/T from T = to 
OO using cr(oo) = a{T*) + J^^, dT (C/T), where we set 
T* = 0.03 Jo . The first term is obtained using Eq. ^) 
and the integral from T* to oo is determined numerically, 
using the QMC results and a high temperature expan- 
sion up to third order. We find a{oo) = 0.68 ± 0.01, in 
very good agreement with the expected result In 2. The 
fraction a{T*) of the entropy is large, contributing ap- 
proximately 25% of the total. Therefore the fact that we 
find the correct value for cr(oo) is a convincing test for the 
correctness of Eq. ( p^ ) , and of our statistical treatment 
of the scaling regime. 

In the study of the random FM-AF spin-1/2 chain, the 
powerful CT QMC loop algorithm allowed us to reach 
temperatures where clear characteristics of the universal 
low-energy properties of this system are observable. The 
characteristic temperature where the scaling behavior be- 
gins depends on the initial distribution of bond strengths. 
Our choice of a non-singular distribution has the advan- 
tage that this temperature is not so low. The key feature 
for analyzing the low-temperature data lies in the fact 
that the spins are correlated within clusters whose length 
grows with decreasing temperature. The scaling of the 
average length / = A^(T/Jo)~^" and the statistics of the 
effective spin distribution allow a very good analysis of 
the numerical data. We determine both the prefactor and 
the exponent. The latter is assumed to be universal and 
agrees surprisingly well with the value found by RSRG 
[0,12|. Among the (unfrustrated) random spin chains 
with complete spin rotation symmetry there are only two 
distinct classes, the ones with the random singlet phase 
fixed point and the others which belong to the class stud- 
ied here [16[ . For the former it was shown by Fisher that 
the RSRG scheme describes the fixed point behavior ex- 
actly pQ]. This was not possible so far for the random 
FM-AF spin chain. Our numerical and statistical analy- 
sis, however, demonstrates very convincingly the validity 
of the universal scaling assumption. Finally, we would 
like to emphasize that our statistical fitting procedure is 
very suitable and useful to analyze not only numerical, 



but also experimental results of this class of random spin 
systems. 
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